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q ■ Abstract 

We study the small-scale behavior of generalized two-dimensional turbulence governed by a fam- 

t— I ■ 

ily of model equations, in which the active scalar 9 = (— A)"/ 2 ^ is advected by the incompressible 
flow u = (—ip y ,ip x ). Here is the stream function, A is the Laplace operator and a is a positive 

i : 

variance (9 2 ) is preferentially transferred to high wave numbers (direct transfer). As this transfer 
O ' proceeds to ever-smaller scales, the gradient V# grows without bound. This growth is due to the 

C0 ■ 

£>y stretching term (V# • V)m whose "effective degree of nonlinearity" differs from one member of the 

Qh. family to another. This degree depends on the relation between the advecting flow u and the active 

I \ scalar 9 (i.e. on a) and is wide ranging, from approximately linear to highly superlinear. Linear 

£T) ' dynamics are realized when Vm is a quantity of no smaller scales than 9, so that it is insensitive to 

o : 

the direct transfer of the variance of 9, which is nearly passively advected. This case corresponds 
! to a > 2, for which the growth of V9 is approximately exponential in time and non-accelerated. 

' For a < 2, superlinear dynamics are realized as the direct transfer of (9 } entails a growth in 

o : 

Vm, thereby enhancing the production of V#. This superlinearity reaches the familiar quadratic 
nonlinearity of three-dimensional turbulence at a = 1 and surpasses that for a < 1. The usual 
vorticity equation (a = 2) is the border line, where Vw and 9 are of the same scale, separating the 
linear and nonlinear regimes of the small-scale dynamics. We discuss these regimes in detail, with 
an emphasis on the locality of the direct transfer. 

PACS numbers: 47.27.-i 
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I. INTRODUCTION 



The production of progressively smaller scales, possibly to be limited by viscous effects 
only, in incompressible fluid flow at high Reynolds numbers is a fundamental problem in fluid 
dynamics. This long-standing problem is of genuine interest for obvious reasons. One is that 
the production of small scales plays a key role in the possible development of singularities 
from smooth initial conditions in the three-dimensional (3D) Euler or Navier-Stokes equa- 
tions that govern the flow. Another reason is that in the presence of a large-scale forcing, 
a persistent production of small scales would be crucial to maintain a spectral energy flux 
(direct energy cascade). The realizability of such a steady and viscosity independent flux is 
central to the Kolmogorov theory of turbulence as this would be required to rid the virtually 
inviscid energy inertial range of the injected energy, thereby making it possible for a statis- 
tical equilibrium to be established. This dynamical scenario is either explicitly or implicitly 
assumed to apply to other fluid systems as well, not just the 3D Navier-Stokes equations. 
For example, in the Kraichnan-Batchelor [1-3] theory of two-dimensional (2D) turbulence, 
the dynamics of the mean-square vorticity (twice the enstrophy) are assumed to be synony- 
mous in many aspects to those of the 3D energy. In particular, the enstrophy injected into 
the system at large scales is hypothesized to cascade to a dissipation range at small scales. 
As another example, the mean-square potential vorticity in the quasi-geostrophic geophys- 
ical flow model is believed to behave in a similar manner [4]. Thus "cascading dynamics" 
have been considered universal among fluid systems. 

The evolution of fluid flow is intrinsically nonlinear because of the quadratic advection 
term, which couples all scales of motion. Apparently, this is an underpinning reason for the 
cascade universality mentioned in the preceding paragraph. However, the "effective degree 
of nonlinearity" of the small-scale dynamics is not always quadratic and differs from one 
system to another. The implication is that the presumed cascades would have fundamental 
differences and would not be universal in a strict sense. For an example of the discrepancy 
in the effective degree of nonlinearity among fluid systems, let us consider the respective 
evolution equations for the 3D vorticity u> and 2D vorticity gradient Vui given by 

d t uj + (u ■ V)w = (w • V)u, V • u = (1) 
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and 



d t Vuj + (u ■ V)Vw = wn x Vw - (Vcj • V)«, V • u = 0, (2) 

where w is the fluid velocity and n is the normal to the fluid domain in 2D. The stretching 
term (u> • V)it for the 3D vorticity u> in Eq. (1) is essentially quadratic in u> because the 
velocity gradient Vu is expected to behave as u> on phenomenological grounds. As a conse- 
quence, an explosive 3D vorticity growth from a smooth initial vorticity field is possible, if 
not inevitable [5, 6]. In contrast, the stretching term (Vw-V)u for the 2D vorticity gradient 
Vcj in Eq. (2) is virtually linear in Vw because Vu is well behaved in the sense that the 
mean square vorticity (u 2 ) = (|Vii| 2 ) is conserved. (Note that the rotation term un x Vc<j 
does not affect the amplitude of Vcj.) As a result, the growth of 2D vorticity gradients 
can possibly be approximately exponential in time only, a relatively mild behavior. Hence, 
one would expect profound differences between the (highly nonlinear) 3D vorticity and the 
(nearly linear) 2D vorticity gradient dynamics. A notable example of these differences is 
that in the inviscid limit, the 2D enstrophy dissipation rate vanishes [7-9], whereas the 3D 
energy dissipation rate presumably remains nonzero. Another example is the discrepancy 
in the dependence on the Reynolds number of the number of degrees of freedom in the two 
cases [10]. 

The effective degree of nonlinearity in the above sense differs not only between 2D and 
3D fluids but also among 2D fluid systems. In this study we investigate this varying degree 
among members of a broad family of generalized models of 2D turbulence, first introduced 
by Pierrehumbert, Held, and Swanson [11]. By so doing we extend several previous studies 
[12-15], aiming to unify our understanding of turbulent transfer in physically realizable 
fluid systems. The family's dynamics are characterized by the material conservation of the 
active scalar 9 = (— A) Q / 2 -?/>, whose variance (6 2 ) is preferentially transferred to high wave 
numbers (small scales). Here ip is the stream function, A is Laplace's operator, and a is 
a positive number. As the transfer of (6 2 ) proceeds to ever-smaller scales, the gradient V9 
grows without bound. This growth is due to the stretching term (V6 ■ V)u whose effective 
degree of nonlinearity depends on a and is wide ranging, from approximately linear to highly 
superlinear. Linear behavior is realized when Vu is a quantity of no smaller scales than 
9, so that the transfer of (9 2 ) to the small scales (direct transfer) has no significant effects 
on Vu. In other words, 6 behaves nearly passively. This case corresponds to a > 2, for 
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which V9 can grow approximately exponentially in time without acceleration. For a < 2, 
superlinear dynamics can be realized as the direct transfer of (9 2 ) entails a growth in Vm, 
thereby enhancing the production of V6>. This superlinearity reaches the familiar quadratic 
nonlinearity of three-dimensional turbulence at a — 1 and exceeds that for a < 1. The 
usual vorticity equation (a = 2) is the border line, where Vu and 9 are of the same scale 
( { | Vxt| 2 ) = (0 2 )), separating the linear and nonlinear regimes of the small-scale dynamics. 
We discuss these dynamical regimes in detail, with an emphasis on the local nature of the 
transfer of (9 2 ). The implication of the present results is that a comprehensive theory for 
this family of generalized 2D turbulence needs to account for the wide range of effective 
degrees of nonlinearity of the family's small-scale dynamics. 

II. GOVERNING EQUATIONS 

The equation governing the evolution of the family of active scalars 9 = (— A) Q//2 -?/> (for 
a > 0) advected by the incompressible flow u = (—ipy,i/j x ) is 

e t + u ■ ve = o. (3) 

This equation has been proposed by Pierrehumbert, Held, and Swanson [11] in an attempt 
to better understand the nature of transfer locality in 2D turbulence, by examining how 
turbulent transfer responses to changes in the parameter a. Equation (3) is physically 
relevant for selected values of a. The usual 2D vorticity equation corresponds to a — 
2. When a — 1, Eq. (3) is known as the surface quasi-geostrophic equation and governs 
the advection of the potential temperature, which is proportional to 9 = (— A) 1 / 2 -*/*, on 
the surface of a quasi-geostrophic fluid. In addition to the genuine interest due to this 
physical significance [12-20], the surface quasi-geostrophic equation has received some special 
attention for its resemblance to the 3D Euler system [21-24]. A mathematical feature of 
particular interest is the possible development of finite-time singularities (from smooth initial 
conditions), which, as argued by pioneering studies [21, 22, 25] of this problem, could be 
associated with the formation of weather fronts in the atmosphere. This, however, appears 
not to be the case [26]. 

For simplicity, we consider Eq. (3) in a doubly periodic domain of size L, and all fields 
concerned are assumed to have zero spatial average. This allows us to express the stream 
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function as 

ip(x, t) = ^2ip(k,t) exp{ik ■ x}. (4) 
fc 

Here fc = 2nL~ 1 (k x , k y ), where k x and k y are integers not simultaneously zero. The reality 
of ip requires ip(k,t) = t/j*(—k,t), where the asterisk denotes the complex conjugate. The 
fractional derivative (—A)"/ 2 (which can be readily extended to a < 0, though not considered 
in this study) is defined by 

6(x, t) = (-A) a/2 ^j{x, t) = J2 ka ^( fc > ex P{^ fc ' x ) 

fc 

= ^9(k,t)exp{ik-x}, (5) 
fc 

where k = \k\ is the wave number. Equation (3) expresses material conservation of 6, which 
gives rise to an infinite set of conserved quantities. In particular, the generalized enstrophy 
(active scalar variance) 

z = \{e 2 ) = \(\{-^) a,2 i>\ 2 ) = \zZ k2a l^( fc > *)l 2 ( 6 ) 

fc 

is conserved. In addition, the generalized energy 

fc 

is also conserved. Note that E is the usual kinetic energy when a — 2, while Z is the usual 
kinetic energy when a = 1. Only for these cases is the kinetic energy conserved. The modal 
powers (spectra) of E and Z differ by the factor k a . Therefore, the redistribution of a non- 
negligible amount of E to small scales would violate the conservation of Z. Similarly, the 
redistribution of a non-negligible amount of Z to large scales would violate the conservation 
of E. This means that if a spectrally localized profile is to spread out in wave number space, 
most of E and Z get transferred to large and small scales, respectively. This is the basis 
for the dual cascade hypothesis in 2D turbulence. Here we are mainly concerned with the 
direct transfer of Z. A more complete treatment should include the inverse transfer of E as 
well since these are known to be intimately related. 

Given Eq. (4), we can express u = (—ip y ,ilj x ) in terms of a Fourier series in the form 

u(x, t) —i S ^(—k y , k x )ip(k, t) exp{ik ■ x}. (8) 
fc 
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By substituting Eqs. (5) and (8) into Eq. (3) we obtain the evolution equation for each 
individual Fourier mode 9{k,t) = k a ^p(k,t) of the conserved quantity 9 

d^,, \ (m a — £ a ) £ x m ^. „ . 

Jt e(k,t)= 0(£,t)9(m,t), (9) 

l+m=k 

where £ x m = £ x m y — £ y m x . The sum on the right-hand side of Eq. (9) involves all 
modes (except 9(k,t)) and is a measure of the level of 'excitation' of the mode 9(k,t) due 
to all admissible wave vector triads k = £ + m. For a given triad, the coupling coefficient 
{m a — £ a ) £ x m/(£ a m a ) depends on a. Its magnitude, together with the magnitudes of 
the coupling coefficients in the governing equations for 9(£,t) and 6(m,t), is a measure 
of triad dynamical activity, in the sense that larger (in magnitude) coupling coefficients 
correspond to more intense modal dynamics. This is intimately related to the effective degree 
of nonlinearity and locality of the small-scale dynamics as will be seen in the subsequent 
sections. 



III. EFFECTIVE DEGREES OF NONLINEARITY OF THE SMALL-SCALE DY- 
NAMICS 

We now examine the behavior of V0. Generally speaking, any derivative (—A) v 9, for 
rj > 0, can be called a small-scale quantity. Here we consider V6 1 , which is a "twin brother" 
of (— A) 1 / 2 ^, for its special status in Eq. (3) as well as its mathematical tractability. For 
a = 2, a similar treatment of A9 = —Au) can be carried out in the same manner. 



A. Growth of the active scalar gradient 

The governing equation for V6* is 

d t V9 + (u ■ V)W = V x uxV9 - (V9 ■ V)u, (10) 

which can be obtained by replacing u in Eq. (2) by 9. Like Eq. (2), the effect of the first 
term on the right-hand side of Eq. (10) is to rotate V# without changing its magnitude. The 
amplification of V6> is due solely to the stretching term (V9 ■ V)w and is governed by 

X79 

d t \V9\ + (u ■ V)|V0| = -T=^ ■ (V0 • V> < |Vu||V0|. (11) 
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Equation (11) implies that following the fluid motion, |V#| can grow exponentially in time 
with an instantaneous rate bounded from above by |Vw|. Hence the behavior of |Vw| holds 
the key to understanding the dynamics of V9. Evidently, following the trajectory of a fluid 
"particle" starting from x = x at t = 0, the growth of |V#| is formally constrained by 

|V0| < |V0 O | exp I jf |V«| drj , (12) 

where 9 = 9(x , 0) and the integral is along the trajectory in question. Hence, on average, 
the rate r defined by 

1 /"* 

r = - \Vu\dr (13) 
t Jo 

provides an upper bound for the exponential growth rate of |V0|. Note that for a — 1 
((|Vw| 2 ) = (|W| 2 )), a double exponential growth of |V0| is allowed but not necessarily 
implied by the preceding equations. Nevertheless, it is interesting to note that Ohkitani and 
Yamada [24] observed such a behavior in their simulations, thereby suggesting a negative 
answer to the question of finite-time singularities in the surface quasi-geostrophic equation. 
This is consistent with the proof of nonexistence of blowup by Cordoba [26]. 



B. Linear versus nonlinear growth of V0 

Now for a sense of the behavior of r, we consider ([Vitl 2 ) 1 / 2 , which bounds (|Vw|) from 
above by the Cauchy-Schwarz inequality (|Vw|) < (\Vu\ 2 ) 1 ^ 2 . For a G [2,4], (\Vu\ 2 ) 1 ^ 2 can 
be estimated in terms of the inviscid invariants using the following version of the Holder 
inequality (see, for example, §5 of Ref. [14]) 

([Vul 2 ) 1 / 2 < (|(-A) a /V| 2 ) 1_2/a (l(-A) Q/2 ^| 2 ) 2/a - 1/2 

= E l - 2/a Z 2/a ' 1/2 . (14) 

So ([V^l 2 ) 1 / 2 is controlled by the inviscid invariants E and Z. For a [2,4], inequality 
(14) reverses direction. Furthermore, if an initial distribution of 9 is to forever spread out 
in wave number space, (|Vu| 2 ) 1//2 increases without bound for this case. This implies that 
there exist different regimes of a for which Vu evolves quite differently, and the active scalar 
gradient dynamics can be characteristically distinct. We discuss all these regimes in what 
follows. 
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For a < 2, the divergence of (IVw) 2 ) 1 / 2 entails an accelerated growth of W from an 
exponential one. This is the superlinear regime discussed in the introductory section. This 
superlinearity reaches the usual quadratic nonlinearity of 3D turbulence at a — 1, where 
(|Vit| 2 ) = (|W| 2 ). Hence, the surface quasi-geostrophic and 3D Euler equations are anal- 
ogous in this aspect. However, the analogy appears to be superficial as the surface quasi- 
geostrophic equation turns out to be far more "manageable" than its 3D counterpart — a 
consequence of the material conservation of 9. For example, a number of global regularity 
results have been proved for the surface quasi-geostrophic equation, by making use of mild 
dissipation mechanisms represented by (— A) n with r] > 1/2 [27-30], which can be much 
weaker than the usual viscosity. Whereas for the 3D Navier-Stokes system, viscosity ap- 
pears to be inadequate for the same purpose. For a < 1, this quadratic nonlinearity is 
surpassed as the ratio (|Vu| 2 )/(|W| 2 ) diverges in the limit (|V0| 2 ) — > oo because 

(|W| 2 ) 2 - Q < (|Vw| 2 )(^ 2 ) 1 ~ a (15) 

(cf. Ref. [14]). Active scalar gradient production can then become highly intense. 

For a G [2, 4], Vu is well behaved in the sense that its mean square is bounded from above 
in terms of the inviscid invariants, [see Eq. (14)]. In this case, Vn is virtually unaffected by 
the direct transfer of (9 2 ). At large t, a general fluid trajectory is likely to have traversed the 
domain many times. The time average in Eq. (13) may therefore be approximately replaced 
by the spatial average. Hence, we can write 

r « (|V«|) < (|Vn| 2 ) 1/2 < £i-2/« Z 2/«-i/2 5 (16) 

where we have used the Cauchy-Schwarz inequality and Eq. (14). This approximation of r 
means that V# can grow exponentially in time without acceleration. Thus, approximately 
linear small-scale dynamics can be expected. Note that 9 behaves almost as a passive scalar 
in this regime. The analogy between this case and that of a passive scalar was suggested by 
Schorghofer [12] on phenomenological grounds. 

When a > 4, inequality (14) reverses direction, and (\Vu\ 2 ) l l 2 can no longer be controlled 
by the inviscid invariants. However, unlike the case a < 2, for which dVu] 2 ) 1 ^ 2 diverges 
toward small scales, when a > 4 velocity gradients can be produced at increasingly large 
scales only. This production depends on the inverse transfer of the generalized energy E 
(Tran 2004). Within the direct transfer range, i.e. the generalized enstrophy range, the 

8 



portion of (|Vtt| 2 ), say Q, cannot increase and instead remains bounded from above in 
terms of Z. More precisely, as the spectra of (|Vw| 2 ) and Z differ by the factor k 2a ~ 4 , we 
have Q < 2kf~ 2a Z (Poincare type inequality), where is the lower wave number end of 
the generalized enstrophy range. This suggests that no significant changes in the effective 
degree of nonlinearity of the small-scale dynamics occur when a exceeds 4. Thus, we can 
expect approximately linear small-scale behavior for all a > 2. 

In passing, it is worth mentioning that while the small-scale dynamics appear to be 
insensitive to a in the regime a > 2, the large-scale dynamics can vary dramatically. The 
reason is that for large a, u is prone to divergence toward large scales as the inverse transfer of 
E proceeds. This undoubtedly intensifies motions at large scales. One may adapt the present 
notion of degree of nonlinearity for a quantitative measure of the large-scale dynamics. 
Analogous to the traditional problem of regularity, which is concerned with the possible 
divergence of W, there is a potential problem that u becomes divergent for sufficiently 
large a if the fluid is unbounded. This interesting problem is left for a future study. 

IV. LOCALITY OF THE SMALL-SCALE DYNAMICS 

This section is concerned with the small-scale dynamics at the modal level. We establish 
a connection between the degree of nonlinearity and dynamical activity of typical local triads 
at small scales. Here the dynamical activity of a given triad is associated with the magnitude 
of the coupling coefficients within the triad and is independent of the amplitude of the three 
modal members. These local triads are shown to be highly active for a < 2 and moderately 
active for a = 2, but become virtually inactive for a > 2. This implies that higher effective 
degrees of nonlinearity correspond to more dynamically intense local triads. Thus, the 
effective degree of nonlinearity is also a measure of dynamical activity of local triads at 
small scales. The transition at a = 2 from high activity to virtually no activity of local 
triads is consistent with phenomenological arguments [11] that the generalized enstrophy 
cascade is spectrally local for a < 2, but becomes dominated by nonlocal interactions for 
a > 2. Below, we also examine the dynamics of nonlocal triads and elaborate on the nature 
of the locality transition, in order to provide a detailed picture of the direct transfer of (9 2 ) 
at the modal level. 

Within each individual triad k = I + m, the transfer of modal generalized enstrophy is 
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governed by 

J t \m\ 2 = Z ££ X m (0(£)e(m)9*(k) + 9*(£)9*(m)9(k)) 
= C k (9(£)9(m)9*(k) + 9* (£)9* (m)9(kf) , 

= C t (9(k)9*(m)9*(£) + 9* (k)9(m)9(£^j , 

j/(™)\ 2 = r ' k p k t Xm (0(k)9*(£)9*(m) + 9* (k)9(£)9(m)) 

= C m (9(k)9*(£)9*(m) + 9* (k)9(£)9(mj) , (17) 

where we have used the identities £xm = £xk = kxm and suppressed the time variable. 
It is well known that both E and Z are conserved for each individual triad. This can be 
readily verified by the fact that the coupling coefficients in Eq. (17), C k , Cg and C m , satisfy 



k a £ a m a 

Furthermore, the transfer of E and Z is from the intermediate wave number to both the larger 
and smaller wave numbers or vice versa (note the signs of the coupling coefficients). The 
former behavior appears to have been observed in numerical simulations of 2D turbulence 
without exception. 

We now analyze the coupling coefficients C k , Cg and C m in detail. As crude estimates 
that hold in general, these can be bounded by (assuming k < I < m) 

= |K-f)txm| k(1 
= r-m°)<xm| 





m a i a 


\(k a 


— m a ) £ x m| 




k a m a 


\(t° 


-k a )£x m\ 



\C m \ = ri [ < (k ■ (18) 

£ a k a 

where we have used \£ x m\ = \£ x k\ < ki. Similar estimates were obtained in [20] 
(for a — 1,2) and in [31] (for a = 1). For a > 2, local triads (i.e. k< £< m) at small 
scales are effectively "turned off" because all C k , C(_ and C m tend to zero in the limit 
k — > oo. Furthermore, the convergence is as rapid as k 2 ~ a . An immediate interpretation 
of this observation is that local triads can be relatively ineffective in the direct transfer of 
(9 2 ) compared with their nonlocal counterparts (see below). At the critical value a = 2, 
C fc , Ci and C m can remain order unity for local triads that satisfy \£ x m\ k 2 and 
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\m a — £ a \ ps \k a — m a \ ps |£ a — fc a | ps A; a . A majority of local triads satisfy both of these 
conditions. They are neither ultra "thin" nor nearly isosceles and correspond to relatively 
sharp estimates in Eq. (18), which reduce to \C k \ ps \Cg\ ps \C m \ pa 1. This means that local 
triads at small scales in the usual vorticity equation are moderately active. They can play a 
significant role in the direct transfer. Finally, for a < 2, the interaction coefficients of these 
same triads diverge as k — > oo. Their divergence can be seen to be as rapid as k 2 ~ a . This 
result suggests that for this case, local triads can play an overwhelmingly dominant role in 
the direct transfer. 

Next, we turn to nonlocal triads. These are "thin" triads with the wave numbers k, £ 
and m satisfying k <C £< m. For this case, C^, Ci and C m can be estimated as follows. 



\Ck\ 


\(m a 


-£ a )£x ml 


a k 2 




m a £ a 




\Ci\ 


\(k a - 


- m a ) £ x ml 


t ik l ~ a 




k Q m a 


\C m \ 




- k a ) I x ml 


ik x - a . 




£ a k a 



(19) 



In the limit t — > oo (while k < oo), vanishes, but both and C m (Ci ~ -C m ) diverge 
as rapidly as £. This implies a vigorous exchange of generalized enstrophy between the two 
neighboring wave numbers £ and m, mediated by a virtually non-participating distant wave 
number k. This ultra local transfer by nonlocal interactions is virtually independent of a 
as the divergence of Ci and C m is insensitive to a. This result implies that local transfer 
by nonlocal interactions is an intrinsic characteristic of this family of 2D turbulence models. 
Note, however, that this transfer can be significant only when the spectrum of the generalized 
enstrophy is not steeper than k^ 1 [32]. In other words, the generalized enstrophy needs to 
be physically present at small scales in order to facilitate such a transfer. This suggests 
that for a > 2 (recall that local triads are dynamically inactive), the generalized enstrophy 
spectra can plausibly scale as k~ l because steeper spectra are unable to support a non- 
negligible direct transfer. This universal scaling has been suggested by Schorghofer [12] and 
Wantanabe and Iwayama [15]. Their justification is that can be considered as a passive 
scalar, a view in accord with the present analysis. 

In passing, it is worth mentioning that the divergence of Ci and C m in nonlocal triads is 
probably the reason for numerical instability in simulations of 2D turbulence with inadequate 
diffusion because local triads with coupling coefficients of order unity are evidently well 
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behaved. Support for this claim can be derived from common observations that numerical 
divergences occur as soon as the modes in the vicinity of the truncation wave number are 
excited and well before they acquire any considerable amount of enstrophy. The same 
instability problem persists for a > 2 although the weak activities of local triads in this case 
may reduce the severity of the instability to a certain extent. 

V. CONCLUDING REMARKS 

We have introduced the notion of effective degree of nonlinearity to quantify the small- 
scale dynamics of a family of generalized models of two-dimensional turbulence governed by a 
broad class of nonlinear transport equations. Here, the active scalar 9 = (—A) a ^ 2 ip (a > 0) is 
advected by the incompressible flow u = (—ip y , ip x ), where is the stream function. We have 
argued that although the advection term is quadratic, the effective degree of nonlinearity of 
the small-scale dynamics is not always quadratic and depends on a. It has been found that 
the active scalar gradient dynamics are virtually linear for a > 2 and become nonlinear for 
a < 2. Furthermore, the degree of nonlinearity increases as a is decreased from 2, becoming 
quadratic at a — 1 and exceeding quadratic nonlinearity for a < 1. It is conceivable that 
credible theories of the family's dynamics, particularly those involving small scales, need to 
account for the dependence on a of the effective degree of nonlinearlity. 

We have also found that local triads at small scales are highly active for a < 2, moderately 
active for a = 2 and virtually inactive for a > 2. On the other hand, nonlocal triads are 
characterized by a vigorous exchange of generalized enstrophy between pairs of neighboring 
wave numbers, mediated by the third non-participating distant wave number. This property 
is common for all a, thereby implying that nonlocal interactions (but ultra local transfer) 
can be considered universal. In the absence of local triad activity (a > 2), this ultra local 
transfer is responsible for the direct transfer of generalized enstrophy. This is similar to 
the problem of passive scalar transport by a large-scale flow as the weak feedback on the 
advecting flow by the active scalar can be neglected [32]. In this case, it appears plausible 
that generalized enstrophy spectra scale as fc -1 . 

The local nature of the generalized enstrophy transfer can be seen to be unambiguous 
in the present study. In general, this transfer is local in wave number space regardless of 
what types of triads make the most contribution. For local triads the generalized enstrophy 
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transfer is inherently local. For nonlocal interactions, the transfer is even "more" local, 
having a relatively higher degree of locality compared to the transfer by local triads. More 
importantly, the transfer between distant wave numbers is largely insignificant. Hence, it 
makes sense to speak of the degree of locality of the direct generalized enstrophy transfer 
rather than to distinguish between local and distant transfer. 
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